Systematics of heavy-ion fusion hindrance at extreme sub-barrier energies 
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The recent discovery of hindrance in heavy-ion induced fusion reactions at extreme sub-barrier 
energies represents a challenge for theoretical models. Previously, it has been shown that in medium- 
heavy systems, the onset of fusion hindrance depends strongly on the "stiffness" of the nuclei in 
the entrance channel. In this work, we explore its dependence on the total mass and the Q-value of 
the fusing systems and find that the fusion hindrance depends in a systematic way on the entrance 
channel properties over a wide range of systems. 

PACS numbers: 24.10.Eq, 25.70.Jj 



During the last thirty years of sub-barrier fusion stud- 
ies, three important observations have been made: 1) 
the discovery of sub-barrier fusion enhancement associ- 
ated with couplings to the intrinsic excitations of the 
participating nuclei 0, 0> EH ; 2) measurements of the 
spin-distributions of the fused-compound nuclei and their 
theoretical description ; 3) the introduction of the con- 
cept of barrier distributions and their subsequent de- 
tailed measurements @, • In these studies it has been 
found that new representations of the fusion cross sec- 
tions, such as the spin distribution da(l)/dl, the mo- 
ments (I) and (I 2 ) of this distribution, as well as the 
quantity d 2 (Ea)JdE 2 associated with the distribution of 
fusion barriers |6(, are essential for exposing pertinent 
features of the data. In general, the coupled-channels 
theory, when using appropriate ion-ion potentials, is able 
to describe the fusion cross section for moderately heavy 
systems down to an energy just below the interaction 
barrier. Recently, it was found that to reproduce fusion 
cross sections at above barrier energies it was necessary to 
increase the diffuseness parameter to values larger than 
those derived from elastic scattering data This ef- 
fect is possibly associated with the opening of the deep- 
inelastic reaction channel. For very heavy systems, it 
is well known that dynamical effects hinder the forma- 
tion of a compound nucleus leading to more complicated 
exit channels such as deep inelastic and quasi-fission pro- 
cesses. The dynamical hindrance of fusion in such sys- 
tems has been described as a diffusion process which may 
eventually reach the configuration of the compound nu- 
cleus 0. 

Recently, a new phenomenon of hindrance in heavy- ion 
fusion reactions has been found in medium-heavy systems 
[ToL ITU If3 , . This hindrance occurs at extreme sub- 
barrier energies whereas the fusion cross section at near 
barrier energies agrees fairly well with standard coupled- 
channels calculations. At present, the exploration of this 
hindrance phenomenon is only in its initial stage; the un- 
derlying physics reason is still unknown. Several colliding 
systems have been measured down to very low cross sec- 
tion levels. In addition, many existing data have been 
reanalyzed in order to uncover systematic trends. Thus, 
it has been found that the nuclear structure of the fus- 
ing nuclei plays a decisive role for the onset energy for 
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FIG. 1: Comparison of the logarithmic derivative and S-factor 
representations of the fusion cross section for the systems, 
58 Ni+ 58 Ni [3 and 64 Ni+ 64 Ni [Q. The dashed curves corre- 
spond to a constant S-factor, whereas the solid curves display 
results of coupled-channels calculations. The L(E) data were 
obtained from a fit to the cross sections at three consecutive 
beam energies. 



the hindrance in medium- heavy systems 0, ^| . In the 
present paper, we study the dependence of the hindrance 
on the mass, and by extension also on the Q-value of the 
fusing systems over a wide range of projectile-target com- 
binations. In general, the Q-value becomes less negative 
with decreasing mass, and even positive for the lightest 
systems. We note that the possible occurance of fusion 
hindrance in the lightest nuclei is of great astrophysical 
interest. 

In order to be able to recognize the hindrance in 
the rapidly varying sub-barrier fusion cross sections, we 
have earlier studied the effect in terms of two represen- 
tations, which are not often used in heavy-ion fusion 
studies. These are the logarithmic derivative, L(E) = 
d\n(aE)/dE and the S-factor, S{E) = oE exp(27n?), 
where r\ = Zi^e 2 / '{Tiv) is the Sommerfeld parameter 
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[ToL ITU 0, 0. In Fig. n these quantities are given for 
the 58 Ni+ 58 Ni and 64 Ni+ 64 Ni systems. The maximum 
in the S- factor occurs at the energy, E s , where the loga- 
rithmic derivative L(E) crosses the curve for a constant 
S-factor, which is given by \H\ 

L^^TTrj/E^OAmZ^^JI/E 3 / 2 (MeV" 1 ), (1) 

where n = AiA2/(Ai+A2) (dashed curves in Fig. [JJ. We 
note that the logarithmic slope of the data, L(E), inter- 
sects L CS (E) at a substantially larger angle, and therefore 
the peak in the S-factor is narrower, for 58 Ni+ 58 Ni than 
for 64 Ni+ 64 Ni. This is a consequence of the "stiffness" 
of the former system. A dependence on the "stiffness" 
of the fusing nuclei has been seen in many cases such as 
90 Zr+ 90 ' 92 Zr, 89 Y (see Table I column 5) and as discussed 
in Ref. ESEj. 

A negative fusion Q-value requires that there be a max- 
imum of S(E) [ll|. This is a consequence of the fact 
that the cross section must vanish at a finite center-of- 
mass energy corresponding to the ground state of the 
fused system, i.e., at E = —Q. In this limit, L(E) = 
a~ 1 da/dE + 1/E — > +oo, whereas L CS (E) remains fi- 
nite. This means that for such systems there is always 
an energy for which the S-factor has a maximum. This 
occurs when L{E) becomes equal to nrj/E, a condition 
that is always fulfilled, since irrj/E is finite near and above 
E = —Q, whereas L(E) — d\n(aE) /dE approaches in- 
finity as E — > —Q. 

For positive Q- value systems, however, a maximum 
may not develop, because both L(E) and L CS (E) become 
infinite in the limit of E = 0. If L(E) does not grow 
faster than L CS {E) with decreasing energy, it may not 
cross L CS (E) for any positive value of E. It is, therefore, 
of interest to study the systematics of the sub-barrier 
fusion hindrance over a wide range of systems, includ- 
ing some with positive Q-value, as is the case mainly in 
fusion between lighter nuclei. 

The expected dependence on the Q-value of the sys- 
tem appears to be borne out by data. The systematics 
of the logarithmic derivative L(E) of fusion excitation 
functions is illustrated in Fig. [2] for a number of systems 
ranging from 10 B+ 10 B to 90 Zr+ 92 Zr. The logarithmic 
derivatives are represented by open circles for five-point 
derivatives, whereas the open squares were obtained by 
a fit to three consecutive data points. We observe that 
L{E) for all systems increases with decreasing energy. 
The dashed curves represent the logarithmic slopes cor- 
responding to a constant S'-factor (Eq. 1). In an earlier 
study of fusion between "stiff" nuclei which did not 
include systems lighter than 16 0+ 144 Sm, we found that 
the S'-factor maximum systematically occurred at a value 
of L s =2.33 MeV -1 corresponding to 



E r /f = 0.356 (ZiZ 2 Vj£)* (MeV). 



(2) 



Studying the full range of systems, we observe that the 
crossing point, E s , for lighter systems, which have in- 
creasingly positive Q-values, indeed occurs at larger val- 
ues of L(E). For the lightest systems, the logarithmic 
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FIG. 2: Logarithmic derivative representations for a range 
of systems. The dashed curves correspond to a constant S'- 
factor, whereas the dashed-dotted curves display results of 
coupled-channels calculations, and the solid curves represent 
a fit to the data using the function a + b/E 3/2 . The range 
of E s values are indicated by vertical line segments for heavy 
systems. For the four lightest systems only lower limits of L s 
can be derived (shown as arrows). The data are taken from 
Refs.: 10 B+ 10 B M, n B+ 12 C , 12 _C_+ 13 C, 12 C+ 16 Q, 16 0_+ 16 



s Ca+ 48 Ca 113, bu Ni+ 8a Y 



and 90 Zr+ 92 Zr 0. 



derivatives of the data intersect the constant S-factor 
curve at a small angle and it is, therefore, difficult to ac- 
curately estimate E s . Consequently, we have used fits to 
the data with the expression a + b/E 3 / 2 (solid curves), 
a and b being adjustable parameters, to obtain a less 
subjective estimate of E s . The results are given in Ta- 
ble I. Relatively large error bars are, however, assigned 
to the resulting E s values and, for the lightest systems, 
only upper limits are given, because of the inaccuracy of 
this procedure. We also observe that the value of the 
logarithmic slope, L(E), obtained by coupled-channels 
calculations for heavy systems (dashed-dotted curves in 
Fig. 2) saturates at a value of ~ 1.5 - 2.0 MeV -1 , much 
lower than measured. It has been shown that coupled- 
channels calculations using reasonable ion-ion potentials 
are unable to reproduce the extreme sub-barrier behav- 
ior [13. 

The systematics of sub-barrier hindrance is illustrated 
in Fig. 3. Here, the derived values of E s and L s = L(E S ) 
are plotted as a function of the parameter Z\ZiyJJL in 
panels a) and b), respectively. Aside from local devia- 
tions of L s from the value of 2.33 MeV -1 in medium- 
heavy systems (of the order of ~ 10%, arising from nu- 
clear structure effects) L s clearly starts deviating from 
this value in lighter systems. The corresponding E s val- 
ues also fall below the E r s ^ systematics (solid curve) 
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TABLE I: This table lists the parameter ZxZ-z^JJl, the energy E s and the logarithmic derivative, L s = L(E S ), that characterize 
the maximum of the S- factor for different systems. Also given are the (dL/dE) exp and (dL/dE) cs , corresponding to the 
measured and the constant 5-factor curves at E s . R is the ratio (dL/dE) exp /(dL/dE) cs , Q is the fusion Q-value, and Vs aS s is 
the height of the Bass barrier H?j . Systems in categories I and II exhibit a clear maximum in the S(E) curve for "stiff' and 
"soft" systems, respectively. A maximum has not quite been reached for systems in category III and IV. Extrapolated values 
of E a and L s etc. are listed for category III, and only upper limits for E s (lower limits of L s ) are included for most category IV 
systems. Uncertainties are given in parentheses for E s , (dL / dE) exp and R. In cases where only upper limits for E s can be given 
the values of (dL/dE) exp , (dL/dE) cs , and R correspond to the crossing points obtained from the fit to the data. Uncertainties 
for L s and (dL / dE) ca can be obtained from the uncertainties on E 3 with the constant S- factor formula. 



System 




E s 
(MeV) 


(Me^ 1 ) 


(dL/dE) exp 
(MeV~ 2 ) 


(dL/dE) cs 
(MeV- 2 ) 


R 


Q 

(MeV) 


(MeV) 


Ref. 


Category I 


yu Zr+ 9U Zr 


10733 


175(1.8) 


2.29 


-1.61(0.16) 


-0.020 


81.9(8.2) 


-157.35 


195.3 


[18] 


90 Zr+ 89 Y 


10436 


171(1.7) 


2.31 


-1.12(0.08) 


-0.020 


55.1(4.4) 


-151.53 


190.1 


ns] 


90 Zr+ 92 Zr 


10792 


171(1.7) 


2.40 


-0.84(0.07) 


-0.021 


39.0(3.6) 


-153.71 


184.4 


[18] 


58 Ni+ 58 Ni 


4222 


94(0.9) 


2.29 


-1.64(0.31) 


-0.036 


44.9(8.6) 


-66.122 


102.0 


[14] 


60 Ni+ 89 Y 


6537 


123(1.2) 


2.38 


-0.80(0.19) 


-0.029 


27.5(6.6) 


-90.497 


136.5 


[10] 


32g_|_89y 


3026 


72.6(0.7) 


2.42 


-0.58(0.15) 


-0.050 


11.5(3.0) 


-36.597 


79.8 


[20] 


Category II 


b4 Ni+ luu Mo 


7343 


121(1.2) 


2.74 


-0.57(0.09) 


-0.034 


17.0(2.7) 


-92.287 


143.3 


ri3] 


64 Ni+ 64 N j 


4435 


87.3(0.9) 


2.69 


-0.35(0.02) 


-0.046 


7.7(0.5) 


-48.783 


98.1 


ri2] 


Category III 


48 Ca+ 48 Ca 


1960 


48.1(0.9) 


2.90 


-0.59(0.03) 


-0.090 


6.5(0.5) 


-2.988 


50.1 


[17] 


28 Si+ 64 Ni 


1729 


47.3(0.9) 


2.57 


-0.70(0.12) 


-0.080 


8.7(1.7) 


-1.783 


50.8 


[21j 


16 0+ 76 Ge 


930.5 


27.6(0.8) 


3.17 


-0.36(0.05) 


-0.172 


2.1(0.2) 


10.506 


32.5 


[22] 


Category IV 


lu O+ i0 O 

12 c+ 16q 


181.0 
125.7 


7.1(0.8) 
<6.2 


4.7(0.7) 
>4.0 


-1.7(0.2) 
-3.0 


-1.0 
-2.2 


1.7(0.2) 
1.4 


16.542 
16.756 


8.2 
6.0 


[16, 23, 24, 25, 26] 
[16, 27, 28, 29] 


12 C+ 14 N 


106.8 


<5.0 


>4.7 


-4.1 


-3.2 


1.3 


15.074 


5.2 


w 


12 C+ 13 C 


89.9 


<4.0 


>5.6 


-3.9 


-2.7 


1.4 


16.318 


4.3 


rie] 


11 B + 12 C 


71.9 


<3.0 


>6.8 


-8.8 


-7.6 


1.2 


18.198 


3.5 


w 


io B+ io B 


55.9 


<1.9 


>10.6 


-26.6 


-24.4 


1.1 


31.144 


2.9 


[15] 



given in Eq. 2. A purely empirical expression 

Ll mp = 2.33 + 400/(Zi Z 2 Vm) (MeV" 1 ) (3) 

(dashed curve in Fig. 3b) is seen to provide a good ap- 
proximation to the experimental data, and it reproduces 
the asymptotic value of 2.33 MeV -1 observed earlier for 
heavy systems with ZxZ^JJi > 2500. The corresponding 
curve for E^ mp obtained from Eqs. 1 and 3, namely 

El m P = (0A$hZ l Z 2 ^Ji/L e s mv ) 2/:i (MeV), (4) 

is seen also to reproduce the experimental values in Fig. 
3a. These two equations thus represent the overall sys- 
tematics for the onset of sub-barrier fusion hindrance. 
This systematics appears to be correlated with the pa- 
rameter Z\Zi^JJi in the simple fashion expressed in Eqs. 
3 and 4, but it should be kept in mind that both the Q- 
value and the fusion (interaction) barrier vary smoothly, 
although not quite monotonically, with this parameter. 
Hence, it is not possible to ascertain whether the ob- 
served physical effect of fusion hindrance is associated 
with cither, or with both of these quantities. 

Figure 3c presents the ratio of the logarithmic slopes, 
R = (dL / dE) exp / (dL / dE cs ) , for the data relative to the 



constant S-factor curve. For ZiZ2y/JI > 2000, this ratio 
is substantially larger than unity which means that there 
is a sharp intersection point between the two curves and, 
consequently, a well defined, narrow maximum in the S- 
factor curve. For Z\Zi^/Ji values below about 2000, the 
slope ratio approaches unity which results in a less well 
defined intercept point. For the lightest systems, it ap- 
pears that the logarithmic slope of the data approaches 
the value for a constant S'-factor and the sub-barrier hin- 
drance may well disappear. 

It should be emphasized that L CS (E) equals the loga- 
rithmic derivative in a point charge, pure Coulomb pen- 
etration model, as long as rj is greater than ~10, which 
is always the case in the energy range studied here. A 
relative slope, i? 3> 1, therefore, implies that the fu- 
sion cross section drops more rapidly than predicted in 
a point charge pure Coulomb interaction model, whereas 
a relative slope near unity indicates that the fusion cross 
section decreases at the predicted rate. 

For orientation, it may also be of interest to relate the 
observed values of E s to the fusion (or interaction) barrier 
and the Q-value of the fusion process, as given in Fig. 3d 
and 3e as a function of the parameter ZyZ^JJi. Since the 
cross section must vanish at an energy of E = —Q (for 
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FIG. 3: (a) Experimental values (symbols) and limits (hor- 
izontal bars) are shown as a function of Z\Zi^fJi for E„ 
(panel a), L s (panel b), logarithmic slope ratio, R — 
(dL/dE) exp /(dL/dE) cs &tE s (panel c), ratio E a /VB ass (panel 
d), and E* = E 3 + Q (panel e), which also shows effective 
particle emission thresholds (histogram). Solid lines in pan- 
els a) and b) correspond to L s =2.33 MeV -1 , whereas dashed 
curves represent empirical trends of the data. Open diamonds 
in panel c) represents the slope ratio at the crossing point ob- 
tained from the fit to the data. 



Q < 0) or E = (for Q > 0) , the larger of these two val- 
ues represents a lower bound for E s . On the other hand, 
one may consider the fusion barrier, taken here from the 
Bass prescription 0] , Vsass , as an upper bound for E s . 
In Fig. 3d, the experimental values of E s do not appear 
to have a simple or fixed relation to Vsass'- the onset of 
sub-barrier hindrance, E Sl occurs at an energy between 
5% and 35% below the fusion barrier. Furthermore, rel- 
atively large fluctuations between systems with similar 
values of the parameter Z\Z 2 ^JJi are present, some of 
which are clearly related to the structure of the fusing 
nuclei [3 El- 

A potential cause of fusion hindrance at sub-barrier en- 
ergies could be the rarefication of final states accessible 



in the fused system, which may be expressed as the ratio 
of total width to the spacing of the states, i.e. T tot /D, 
in the appropriate energy regime of the compound nu- 
cleus. In other words, a fusion reaction can only proceed 
if a quantum state with the appropriate energy, spin and 
parity is available in the compound nucleus. The T tot /D 
value is expected to increase exponentially with excita- 
tion energy and will approach unity slightly above the 
particle (n, p, or a) emission threshold (binding energy 
+ Coulomb barrier). In order to explore this possibility, 
we have plotted in Fig. 3e, the excitation energy corre- 
sponding to E s , i.e., E s + Q (solid symbols) and com- 
pared it to the particle thresholds for the systems listed 
in Table I. We note that the experimental values of E s 
correspond to excitation energies exceeding the particle 
thresholds by a wide margin in all cases, and there does 
not appear to be a significant correlation between these 
two quantities. Even accounting for the fact that some 
of the excitation energy is bound in rotational energy 
does not alter this conclusion. Although the simple ex- 
planation of rarefication of the final states in the fusion 
process is appealing, it does not appear to account for 
the observed hindrance phenomenon. It seems that the 
behavior seen in Fig. 3a-c provides indications that the 
hindrance phenomenon is closely related to the entrance 
channel. 

In conclusion, the systematics of sub-barrier fusion 
hindrance has been studied over a wide range of sys- 
tems from 10 B+ 10 B to 90 Zr+ 90 Zr. Hindrance appears 
to be a general phenomenon, at least for systems with 
ZiZ 2 ^/Ji > 3000. For the lightest systems {Z 1 Z 2y /JI < 
200), the logarithmic slopes of the cross section in the 
sub-barrier region merge smoothly into those expected 
on the basis of a constant S'-factor (i.e. a point charge 
in a pure Coulomb interaction model). Simple empiri- 
cal formulae are given for both the energy and the log- 
arithmic slope of the cross section at which the onset 
of fusion hindrance occurs. These point to an entrance 
channel effect as the source of this phenomenon. Un- 
til now, all direct observations of fusion hindrance have 
been made in systems with Z\Zi^fJi > 3000. It would be 
interesting to study sub-barrier fusion in more systems 
in the range 200 < Z\Zi^JJi < 3000, where the fusion 
Q- values change from positive to negative values. If the 
fusion hindrance does indeed occur in light systems, such 
as 12 C+ 12 C, 12 C+ 16 and 16 0+ 16 0, it will strongly af- 
fect the predicted rates of astrophysical processes, which 
are presently obtained by simple empirical extrapolations 
from experimental data. As yet, no satisfactory theoreti- 
cal explanation for this phenomenon has been put forth. 
Simple considerations in terms of relations to the fusion 
barrier height or the rarefication of compound states in 
the fusion channel do not appear to clarify the situation. 

This work was supported by the U.S. Department of 
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W-31-109-ENG-38. 
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